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Abstract 

We study the behavior of the entanglement entropy in (2 + l)-dimensional strongly coupled 
theories via the AdS/CFT correspondence. We consider theories at a finite charge density 
with a magnetic field, with their holographic dual being Einstein-Maxwell-Dilaton theory in 
four dimensional anti-de Sitter gravity. Restricting to black hole and electron star solutions at 
zero temperature in the presence of a background magnetic field, we compute their holographic 
entanglement entropy using the Ryu-Takayanagi prescription for both strip and disk geometries. 
In the case of the electric or magnetic zero temperature black holes, we are able to confirm 
that the entanglement entropy is invariant under electric-magnetic duality. In the case of the 
electron star with a finite magnetic field, for the strip geometry, we find a discontinuity in the 
first derivative of the entanglement entropy as the strip width is increased. 
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1 Introduction 


The AdS/CFT correspondence [1, [2, [3j, 2] has become a prevalent new theoretical tool for un¬ 
derstanding strongly-coupled, non-gravitational physics. Its influence has extended beyond the 
realm of strict high energy physics, like understanding the quark-gluon plasma , into con¬ 
densed matter systems. Of interest are strongly-coupled charged fermionic systems, where it is 
hoped the low-energy physics can be studied via the correspondence (see refs. 0 El El ED]), in 
particular, progress has been made in understanding Fermi Liquids (FL), non-Fermi Liquids (NFL) 
and Fractionalized Fermi Liquids (FL*), phases of compressible metallic states of quantum matter 



These different phases with finite charge density can be characterized in terms of their 
violation/agreement with the Luttinger relation, which relates the total charge density Q to the 
volumes enclosed by the Fermi surfaces at zero temperature. Compressible states of matter with 
finite charge density dual to charged black holes PS ng [m fl8| violate the Luttinger relation. 
However, holographic duals to metallic states, such as the electron star solutions of refs. [ 131221121 ], 
do satisfy the Luttinger relation [22 , 23]. The essential idea [12] is that when the field theory charge 
density Q is dual to gauge-invariant fermions - “mesinos” - in the gravity bulk, their Fermi surfaces 
do satisfy the Luttinger relation and are said to have visible Fermi surfaces. When the charge density 
is instead sourced by a charged horizon, then the Fermi surfaces are said to be hidden , leading to 
a violation of the Luttinger relation and thus describe “fractionalized” charged degrees of freedom. 
In the case when there are both gauge-invariant fermions and a charged horizon in the bulk, like 
the solutions in ref. [21], then the phase is considered to be “partially fractionalized.” 

Among the interesting questions posed about NFL and FL* phases is the issue of how to 
characterize the presence of Fermi surfaces when the conventional methods from field theory are not 
easily applicable. It was shown in ref. El that such phases in the field theory have a logarithmic 
violation of the area law for the entanglement entropy, and this sets a strict criterion for gravitational 
models to be considered dual to such phases. The prescription of Ryu and Takayanagi [25] for the 
calculation of the holographic entanglement entropy allows us to directly address this question. 

Results of holographic computations suggest that the mesonic phases are dual to FLs, 
whereas the fractionalized phases are dual to NFLs and partially fractionalized to FL*s [23]. In 
particular, it was found in refs. nn ei that the hidden, gauge-charged Fermi surfaces in the bulk 
do lead to a logarithmic violation of the area law when the metric is of hyperscaling-violating form 
with particular values of the dynamical critical exponent and the hyperscaling violating exponent. 
In the case of the electron star, the logarithmic violation depends only on the charge sourced by 
the horizon and not that coming from the star itself m- 

Here we extend existing work in the literature by computing the entanglement entropy of 

1 See refs. 11311141 for explanations of FL, NFL and FL* phases, where we have used their abbreviations for the 
names of the phases. 
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solutions to a (3+l)-dimensional Einstein-Maxwell-Dilaton (EMD) theory in the bulk that are not 
of hyperscaling-violating type and additionally have an external, constant background magnetic 
field turned on, including the magnetic electron star solutions found in ref. [26]. 

Magnetic fields are also important probes of the physics of transport (especially for (2 + 1)— 
dimensional systems, to which our gravity solutions are dual), and our work is intended to be 
part of a program of trying to understand, using such external fields, possible new ways to classify 
characteristic behaviours in the phases we capture. In this paper we will focus on studies of the 
behaviour of the entanglement entropy in the presence of the magnetic field. Since, as stated above 
(see also further discussion below) our geometries are not of hyperscaling-violating type, we do 
not expect the characteristic logarithmic violations of the area law, and so we seek to explore and 
exhibit the entanglement entropy’s behavior in the various regimes to which we have access. 

Our work studies gravitational backgrounds not previously studied in the literature; in the 
notation of ref. [12] . their a and /?, which determine the dilaton’s potential and coupling to the 
Maxwell sector, in our cases are equal: a = f3. In terms of the dynamical critical exponent z, this 
corresponds to z —> oo and a hyperscaling violating exponent 9 = 2. These values are incompatible 
with the form of the hyperscaling-violating metric, and in particular do not satisfy the requirements 
of ref. [12] that (5 < and 9 = 1. 

In ref. EH, the authors consider an EMD theory using a potential for the dilaton of the form 
E(<h) = — |Vo|exp (25$) and coupling to the Maxwell sector Z(<f>) = exp(2a<I>). They study the 
entanglement entropy of purely electric, purely magnetic (argued for via electromagnetic duality), 
and dyonic solutions of this system, adding the magnetic field as a small perturbation, B <C /i 2 , in 
their dyonic system. (Here, fj. is the chemical potential.) In their notation, our solutions correspond 
to the choice a = 5, which is a case they do not consider in their woi'10 In particular, they start 
with a hyperscaling-violating metric for the purely electric solution parameterized by their (a, 5), 
but in the case a = 5, the dynamical exponent and the hyperscaling violating exponent both 
diverge, positively and negatively, respectively. Thus, our EMD solutions lie along a line in their 
phase diagram for which they do not explore the entanglement entropy. In addition, our dyonic 
solutions have Q = B. 

We present our work as follows. In Section [2] we begin with a review of the holographic 
entanglement entropy, but present the detailed derivations of the formulae we display in the Ap¬ 
pendix. In Section [3] we present our gravity backgrounds, although for more information on their 
derivation and the magnetic electron star solutions we refer to ref. [26]. Section [4] shows the results 
for the strip and disk (see Section [ 2 ] for a definition of these) entanglement entropy for the purely 
electric dilaton black holes, and Section [5] for the purely magnetic dilaton black holes. Section [6] 
studies the behavior of the strip entanglement entropy under electromagnetic duality for the purely 
electric and purely magnetic dilaton black holes. In particular, it is shown how the entanglement 

2 In ref. [38], the authors consider a fermionic two-point function for a = 5 and find that it exhibits non-Fermi 
liquid behavior. 
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entropy is invariant as long as one is careful to take into consideration the relative positions of the 
physical horizons. The entanglement entropy for the strip and disk for the dilaton-dyon black hole 
is presented in Section [7J and then in Section [ 8 ] we consider the mesonic phase of the magnetic 
electron star solutions of ref. {261 . We end in Section [9] with our conclusions and a discussion of 
future worl 0 


2 Review of the Holographic Entanglement Entropy 


The holographic entanglement entropy of Ryu and Takayanagi [25| computes the entropy of en¬ 
tanglement Sa between two subsystems A and B = A (the complement of A). The prescription 
involves finding the minimal area surface that extends into the AdS gravitational bulk, whose 
boundary at conformal infinity is that of subsystem A. That is, if 7^4 is the minimal surface in the 
{d + 2)-dimensional bulk such that d'jA = A at the (d + l)-dimensional UV boundary, then 


£4 = 


Area(yA) 

(d+2) 


4 G 


( 1 ) 


N 


where the surface 7^4 has co-dimension 2 and Gy 1-2 is Newton’s constant in (d + 2{-dimensions. 
To compute the area, ones takes the bulk metric and integrates its pull-back H^ v onto 7,4 for 
a constant time slice t = fo : 


Area(7,4) = / d d xVJl . ( 2 ) 

Jy a 

Here x represents the coordinates on 7 a, which are generally given as the embedding of the surface, 
and H is the determinant of H^ u . 

The above prescription was described for (d + 2)-dimensional AdS spacetime, but it applies 
more generally. In particular, there may be a non-zero and dynamical dilaton in the higher 
ten-dimensional string theory. When this is the case, eqn. ([ 2 ]) still holds so long as the metric is 
written in Einstein frame IHl- 

In what follows, we will work with a four-dimensional bulk (so that d = 2 in eqn. ({Tj) and 
hence will drop any sub- or superscripts denoting the dimension) with Newton’s constant in eqn. 0 
given by dimensional reduction in the usual way, and the regions we consider in the boundary are 
that of the strip and the disk as seen in Figure [lj we shall denote their entanglement entropies 
as Ss and Sd, respectively. In particular, we will be interested in the finite part of the entanglement 

3 While this work was in preparation, two papers appeared on the arXiv, refs. [291 l30l that construct magnetic 
electron stars that are different from the ones presented in ref. [26]. In this new work, there is no dilaton present 
in the theory and their horizons have finite temperature. They do not consider the entanglement entropy of their 
solutions, but it would be interesting to compute the entanglement entropy for their solutions and compare them to 
the results in this paper. 
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entropy, ss and so, defined by 


4 G n S s 


2 L 2 ££ 



AGjySo = 2 ttL 2 


(^sd + 


(3) 

(4) 


Here «£? denotes the strip geometry’s “infinite” side length, L is the AdS length scale, e is a 
chosen UV cutoff, and £ is the radius of the disk. The added term on the right hand side of 
eqns. ^ and Q removes the leading divergence of the entanglement entropy, and we note that 
these expressions have been defined for dimensionless coordinates so that «£?, e, £, and the ss/d are 
all dimensionless. 
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Figure 1: A schematic representation of the two regions, A (with complement B), on the boundary 
(with coordinates {x,y}) whose holographic entanglement entropy we consider, where the holo¬ 
graphic direction is z. The minimal surfaces hang down into the bulk, (a) The strip, whose finite 
entanglement entropy we denote ss, with finite width £. (b) The disk, whose finite entanglement 
entropy we denote so, with radius l. 

For reference, we now present generic formulae for the finite entanglement entropy for the 
strip and the disk for a form of the bulk metric that all of our geometries will have: 

ds 2 = -g u (z)dt 2 + g zz (z)dz 2 + g xx (z)dx 2 , (5) 

where we have taken the UV to be at z = 0 and the horizon to be at some z = zh that we rescale 
to be at z = 1. Writing the M 2 part of the metric as dx 2 + dy 2 and having the finite length of the 
strip run from — 1/2 < x < £/2, the minimal surface will be symmetric about y = 0. We denote 
the ^ value of the turning point of the surface as zt- Let the surface have coordinates (x,y), the 
same as the M 2 coordinates, with embedding 2 = z(x). 

In the Appendix we derive the formulae for the entanglement entropy for the strip and 
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the disk geometries. We note that in the case of the strip, the second order problem of finding 
the minimal surface z(x) can be reduced to a first order problem via an integral of motion (see 
the Appendix for details), which allows the integrals to be written in the form shown below; in 
addition, it grants more control on the numerics of our work in the following Sections. For the 
strip, the result i^] 

4 G n S s = 2^ [ T dz 9xx{z) 2 == ; ( 6 ) 

£ V Sim ( 9xx- 9xx(z t ) 2 ) 


where we have introduced the UV cutoff e in the lower limit of the integral. The length of the strip, 
1=2 Jq 2 dx, as a function of the turning point is given by 


e 

2 



dz 


9xx {.zt ) 

\! Sim ( 9xx ^z ) 2 - 9xx(z t ) 2 ) 


(7) 


In the case of the disk geometry, the full second order problem must be solved, so the 
formula for the entanglement entropy is simply given by the pull-back 

40kSd = 2 ”l drrs ^ z{T)) f^m (i) • <8) 

where 0 < r < i is the radial variable for polar coordinates of the M 2 , i.e., ds^ 2 = dr 2 + r 2 d0 2 . In 
this case, we must explicitly solve for the minimal surface z(r) and then input that into eqn. (|8]) 
to find the entanglement entropy. Because of this more numerically intensive approach, we will 
observe some increased numerical variance for the disk results relative to the strip in the Sections 
that follow. 


3 Gravity Backgrounds 


Our main goal will be to study the entanglement entropy for the magnetic electron star solution 
found in ref. [26j, but to understand that result in a larger context, we will consider other related 
backgrounds as well. Here we present the common aspects of these solutions, giving the precise 
forms of the metrics below when we compute their entanglement entropy. 

We will begin by considering zero temperature black hole solutions of a four-dimensional 
Einstein-Maxwell-Dilaton system with action 


Semd = J d A x\/—G 


^2 ( R - 2d^d^ - H(4>)) - 


(9) 


4 Note that here we have not yet subtracted the leading AdS divergence nor pulled out a factor of the AdS radius L, 
so that these are not quite the expressions for ss and So- 
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Here k 2 = 8ttGn is the gravitational coupling and e is the Maxwell coupling. We take the potential 
for the dilaton and its coupling to the Maxwell field to be 

V($) = --^cosh(2$/%/3), Z(4>) = e 2 ^ . (10) 

We will consider solutions that have either just electric or just magnetic charge, as well as 
solutions with both. Our ansatz for the metric and Maxwell fields with a constant magnetic field B 
turned on is 

T 2 

ds 2 = — 5 - (-F(z)dt 2 + G(z)dz 2 + A(z)dx 2 ) , (11) 

z z 

F zt = —h'(z ), F xy = . (12) 

KZh KZj] 

The horizon is at zh and the UV at z = 0. When we wish to consider cases where there is no 
electric or no magnetic field, we will find it consistent to simply set h(z) = 0 or B = 0, respectively. 
The coordinates are dimensionful; however, the fields h and B are dimensionless. We will require 
our solutions to asymptote to AdS 4 in the UV, which dictates the leading behavior of the fields to 
be 


F(z) = 1 + ..., G(z) 

h(z) = p — —z + ... , 
Zh 


1 + . 

• • ! — 1 + . . . , 

(13a) 


01 . 02 2 . 

— —z H— ttZ + ... . 

(13b) 

ZH Zfj 


The physical quantities in the dual field theory are the chemical potential pp = p^Pi the charge 
density Qp = -pp-Q, the magnetic field Bp = -^4 B, the source for the operator dual to the 

jj kz h 

dilaton 0 i p = ^|^ 0 i, and the vacuum expectation value (vev) of the operator dual to the dilaton 

02 p = — 4 t 02 - It is convenient to work with dimensionless quantities, and so, by taking ratios with 

kz h 

respect to //p, we can characterize the dual field theory by the following dimensionless ratios: 


Bp k B 
p? p eL n 2 


Qp _ n Q 
Up eL fi 2 


01 P _ 1 01 

HP e 11 


02 P _ K 02 
fi 2 p e 2 L p? ’ 


(14) 


In order not to specify values for e, L, and k, we will give the values of B jp 2 , Q/p?, etc. in what 
follows. We also work in dimensionless coordinates defined by replacing ( 2 , t, x) in the expressions 
above by (zpz, zpt, zpx), such that the horizon position is at z = 1 . 

We now summarize our work in ref. [ 25] on how to introduce the charged star, and refer to 
that reference for the full details. We take as our action 


S = SeMD + Sfluid, 


(15) 
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where the fluid action comes from the Lagrangian IMEIIMI 


Aiuid = 'f I G (~p(a) + au * ( d ^ + A^ + ad^/3) + A (u^u^ + 1)). (16) 


Here (f> is a Clebsch potential variable, (a, (3) are potential variables, and A is a Lagrange multiplier. 
The energy density of the fluid, p, and the charge density, a , were found in ref. |2B] to be that of a 
free fermion with mass m = nm/e: 


1 


a = eL 2 Kcr = (p, — rh 2 ) 

O 


2 n3/2 


p = L 2 K 2 p = ^p,y/p 2 — rh 2 (2p 2 — m 2 ) + m 4 In 



(17) 

(18) 


where p(z) = zh(z) / y/F(z), /3 = e A L 2 (3/n 2 ~ 0(1) is a constant of proportionality, and a tilde 
means that it is a dimensionless quantity. The free fermion result is valid so long as 2 qB <C 1, 
where q = and a(z) = A(z)/z 2 . To be in a regime where we can use classical gravity we 

must have k/L <C 1, and since /? ~ 0(1) it follows that e 2 ~ k/ 1 < 1, Thus the free fermion result 
is valid when <C 1. In the solutions that follow it is assumed that we adjust the value of e so 
as to maintain this constraint. Our ansatz for the fluid is 


eL 


Ut 


= Ay = — By, 


= —Bxy , a = —— By , 8 = x. (19) 

K K 


With this ansatz the equations of motion from the action in eqn. (15) become, after some 


work, 


P\ z ) + (p(z) + p( z )) ~ °( z ) 


2f( z 


h'(z) 

y/Kz) 


= 0 


a"( Z ) 

A*) 

f(z) 


- a (z) ( 


g \ z ) , m 
2 g(z) 2 f(z) 

(g\ z ) , m 

\2g(z) 2 f(z) 


+ a (^) (aiz) (p(z) + p( z i) + 2 $'(-) 2 ) 

2a'(z) A 

a( z ) ) 


= 0 , 


<$>'(z) 2 + g(z) 


+ - 2 &(z) 2 - g(z) ( 5 P(z) + p(z) - 2H(d>))) = 0 , 

z(m 2 , ( P(A 1 _ a 'W( z ) _ z(^)h'( Z ) 2 

2a(z) 2 V 2 ’)) 4 a(z) 2 2 a(z)f(z) 2 f(z) 


(20a) 

(20b) 


(20c) 
0 , (20d) 


h"( Z ) 


h'(z) 


g'jz) /'(*) 

2 g(z) 2 f(z) 


a/(z) 

a(z) 


<P"(z) + <l>'(z) 



g\ z ) , a \ z ) \ 

2 g(z) a(z) ) 


z'(<$>)&(z) \ Vmgjz) 

Z($) ) Z($) ' 

_ Z^} f g(z)B 2 _ h'(z) 2 \ g(z)V'($) 

4 V a( z ) 2 f( z ) J 4 


(20e) 

(20f) 


In these equations, we have written f(z) = F(z)/z 2 , g(z) = G(z)/z 2 , a(z) = A(z)/z 2 , and 




































V(<h) = L 2 V(Q). Here the fluid pressure P(z) is given by the thermodynamic relation P(z) = 
—p(z) + jj,(z)a(z) so that eqn. ( 20 a) is automatically satisfied. As a check, setting the fluid fields to 
zero reduces eqn. (20) to those coming from just the Einsten-Maxwell-Dilaton system in eqn. ([ 9 ]), 
so when we consider cases without the star we will find solutions to eqn. ( 20 ) with the fluid fields 
turned off. We still require our star backgrounds to be asymptotically AdS 4 , and so the UV dual 


field theory is again characterized by the dimensionless ratios in eqn. (14) 
It was found in ref. 


that the solutions of these equations fall into three broad phases 
governed by the IR behavior of the dilaton. In the case with no star, the dilaton can diverge 
logarithmically either positively or negatively in the IR, giving rise to a purely electric horizon 
or a purely magnetic horizon. If it tends to a finite value, then it is a dyonic solution. Figure [2] 
schematically illustrates this phase structure. With the star present, the same classification scheme 
can be used and we can label the solutions: a “mesonic” phase, where all the electric charge is 
sourced by the star; a “partially fractionalized” phase, where a fraction of the charge is sourced 
by the star and the rest by the horizon; and a “fully fractionalized” phase, where all the charge is 
sourced by the horizon. 



Figure 2: Schematic phase diagram of the various black hole solutions considered in this paper. The 
vertical axis is the purely electric dilaton black hole (PED) with no magnetic charge, the horizontal 
axis is the purely magnetic dilaton black hole (PMD) with no electric charge, and the line is the 
dilaton-dyon black hole, with magnetic and electric charge satisfying Q = B. The IR behavior of 
the dilaton is indicated for each phase, in coordinates where the horizon is located at infinity. 

In this paper, we will study the entanglement entropy of the following phases: no star purely 
electric horizon, no star purely magnetic horizon, no star dilaton-dyon, and star mesonic phase. 
These solutions were found in ref. [26) as asymptotic expansions in the IR and were integrated 
numerically out to the UV, matching onto AdS. 4 . We display the relevant expansions as they are 
needed in the following sections. 
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4 Entanglement Entropy of the Purely Electric Dilaton Black 
Hole 


We first begin with the purely electric dilaton black hole (PED) with no star. An analytic form of 
this solution was written down in ref. [35] (where it arises as the three-equal-charge dilatonic black 
hole in four dimensions) at a single point in the phase diagram; below, we indicate which point this 
is in terms of parameters. Since our star backgrounds are constructed numerically, as a check of 
our numerical procedure, we construct a numerical PED below and compare the results with the 
analytic solution of ref. {351 . 

Our IR (z = 1) series expansion for the metric 0 is 


F(z) 

G{z) 

A{z) 


(1 - ^) 3/2 /o 


OO 


l + X^Mi-zy 

n= 1 


( oo 

X>n(l-*) n + <S«7(l 

n=0 

(1 - z ) 1 / 2 . 


+ <5/(1- 




(21a) 

(21b) 

(21c) 


We can consistently set B = 0 in our equations of motion and ansatz, and take for the electric 
component of the Maxwell field and for the dilaton the expansions 


/ *T 

1 

1—1 

'-o 

+ 

1 

r—1 

II 

TT 

(22a) 

/o 

4>(^) = —— log (1 — z) + ^(l — z) b . 

(22b) 


We have turned on a perturbation {5f,5g,5h,5&}, with b = | (—3 + V57 ), to allow us to flow to 
different values of the dual field theory parameters in eqns. ( |14[ ). To find the perturbation, we treat 
the equations of motion (EOM) as functions of the 5f, etc., and then consider the linear problem 


<9EOM(<5/) 


ddj 


Sj=0 


<5j = 0 , 


(23) 


where I, J = 1,...,4 and <5i = 5f, 62 = 5g, etc., and “EOM” in eqn. (23) refers to each of 


the eqns. (20). We require this to have a nontrivial solution for the <5/ and thus, the matrix of 


derivatives must be non-invertable. Requiring its determinant to vanish gives us a polynomial 


condition for /3 which we can solve, and then using this, we can solve the linear system in eqn. (23) 
for the perturbations. 


Substitution of eqn. (21) into the equations of motion in eqt. (20) fixes all coefficients in 


the expansions except for /o, which is chosen so that the UV metric flows to AdS. In fact, this fixes 
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all perturbations except 5f , whose choice takes us to different values of the dimensionless ratios in 


eqn. (14). 

We now compute the finite part of the entanglement entropy for the strip geometry using 
eqns. (J6]) and ([7]). Figure [3] shows a plot of ss as a function of the strip width t. In order to check 
our numerics, we compared this calculation with that of the analytical solution for the PED in 
ref. |35r[ Their metric, written in the form of eqn. (11), is 


F(z') 

G(z') 

A(z') 


1 + Qz 1 


3/2 


/ 


,3 J3 


1 - 


Q A Z 


1 + Qz 1 


A "3/2 


^ ^1 + Qz ' 

( 


3~/3 


1 - 


Q A Z 


^ ^1 + Qz ' 


1 + Qz 1 


A 3/2 


(24a) 

(24b) 

(24c) 


where we have used a prime to denote their coordinates, and in their coordinates {f', 2 :',x'} the 
horizon is at z' = oo and the UV at z' = 0, and L' is the AdS radius. Lastly, to compare to our 
solution we did one further rescaling of the coordinates by Q so that eqns. (24) are only in terms 
of z ', although in Section^ we need the explicit Q dependence and so have kept it in for eqns. (24). 

The analytic solution corresponds to our 5f = 0, for which the dimensionless UV field 
theory parameters take the values 


u = 1.2247 , — = 0.3535 , H, = 0.8165 , % = 0.1444 . (25) 

/r /V /V 

These values agree with the analytical values to within 10~ 4 accuracy. The result of the analytical 
entanglement entropy computation is also shown in Figure [3j and there seems to be a mismatch. 
However, this is due to the difference in coordinates used in the two calculations. The analytic 
solution’s coordinates are related to our coordinates via z = jq-v- The consequence is that when 
regulating the entanglement entropy so that it is finite, we must specify what our chosen regulator 
is and this is a coordinate dependent statement. When comparing the two results we must take into 
account that the UV regulators are different because the location of the UV relative to the IR is 
rather different in the two cases; indeed, the shift seen in Figure [3] for the asymptotic IR value of ss 
is 1, which from the coordinate relationship between the two systems is precisely the difference in 
regulators: l/e = 1 + 1/e'. 

5 In comparing their four-dimensional Lagrangian with ours, the relationships between their dilaton a and ours 
is a = In addition, their metric is in terms of a coordinate r'. but we have done a change of coordinates to 

z' = L' / r' and replaced their {t', x'} by L'{t ', x'} for comparison. We have also defined a dimensionless charge Q via 
Q = Q’/L'. We have omitted writing expressions for their dilaton and Maxwell field as their particular form is not 
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ss 



Figure 3: The finite part of the entanglement entropy ss as a function of the strip width £ for the 
numerical PED (blue/lower), with Sf = 0, and the analytical solution (red/upper), eqns. ([24]). We 
have taken our UV cutoff to be e = ICE 5 and our maximal IR value of zt = 0.9999. See text for 
explanation of the shift. 


It is usually expected that the large £ limit of the entanglement entropy should approach 
the thermal entropy of the system. However, in this case, our PED solutions have zero thermal 
entropy and we see that both the numerical and the analytic solution asymptote to two different 


values as £ increases. In fact, we see in Figure 4(a) plots of the strip entanglement entropy for 
various values of 5f each of which asymptotes to a different large £ value. We also find that as Sf 
is increased, the value of the dual held theory chemical potential g also increases, while the values 
of (f> \/n and Q/fi 2 decrease. 

Interestingly, if we consider the “in-falling” minimal surface, namely the one that just hangs 
straight down into the black hole, and compute its finite area via 


Al ^ /lf ^ = [ dzy/a(z)g(z) - * , 


(26) 


(where y;f denotes the in-falling surface) we find that this value agrees to within reasonable nu¬ 
merical precision with the large £ value of ss- For example, when Sf = 0, we find that the large £ 
value of ss is given by —0.4325 and the finite area of the in-falling piece is given by —0.4427. 
This suggests that these “side pieces” contribute to the entanglement entropy. Since our solutions 
have zero temperature and zero thermal entropy, as well as a zero area horizon, this may be the 
reason we are able to see their contribution, whereas in the case of a non-zero thermal entropy, 
these pieces are negligible compared to the thermal entropy. Figure |4(b)| shows the same plots as 


Figure 4(a) but zoomed in to show a comparison between the large £ behavior of the entanglement 


needed for what we discuss here. 
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(a) 


(b) 


Figure 4: (a) The finite strip entanglement entropy ss as a function of £ for various values of Sf. 
Blue/lower is Sf = —10 ( Q //i 2 = 1.0209), red/middle is Sf = 0 ( Q/n' 2 = 0.8165), and orange/upper 
is Sf = 10 (Q/n 2 = 0.5863). (b) A zoomed in version of (a) where the horizontal lines are the areas 
of the in-falling solution for each. 


entropy and the in-falling piece. In Figure [5j we show how the large £ value of ss changes as the 
dimensionless ratio Q/fJ. 2 is varied. 


ss 
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Figure 5: The large £ value of ss as a function of Q/fi 2 for —10 < 6 f < 10. 

We can also compute the entanglement entropy for the disk geometry s d defined in eqn. Q 
and found by computing the integral in Q using the minimal surface z(r) found by extremizing 
the area functional of ([8]). For numerical reasons, we find it more convenient to work with a new 
radial coordinate defined by [36] 

((£) = i v /£2 + e 2_ r 2 ) ( 27 ) 

where £ is the radius of the disk. In terms of this radial coordinate the pure AdS 4 minimal surface 
would just be given by z(C) = £(. Using this, we are able to find sd for various values of Sf. In 
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Figure [6j we plot sd/£ vs. £ for three different values of 5f. We se^jthat entanglement entropies 
seem to be asymptoting to different large £ values as in the case with the strip. 


sd 

l 



Figure 6: The finite part of the entanglement entropy spasa function of the disk radius £ for three 
values of 5f. Blue/lower is 5f = —4 ( Q/fj? = 0.9002), red/middle is 5f = 0 ( Q //t 2 = 0.8165), and 
orange/upper is 5f = 4 (Q//u 2 = 0.7289). Recall that finding sd requires solving a second order 
differential equation as well as numerical integration, so the numerical accuracy is not as controlled 
as in the strip case, hence the numerical variance in the entanglement entropy curves shown. 


5 Entanglement Entropy of the Purely Magnetic Dilatonic Black 
Hole 


In this section, we will consider dilatonic black holes that have no electric charge but do have a 
magnetic charge. We can use as our ansatz the same metric © and Maxwell fields (|12[) except 


6 As mentioned in Section[2]finding sd requires solving a order differential equation as well as numerical integration. 
Since our backgrounds are numerical and we also have to solve the second order differential equation numerically, we 
have some numerical variance, as seen in the disk plots. 
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we set h(z) = 0. The IR expansion for our metric fields and dilaton then becomes 


/(*) 
9(z) 
a (z ) 

$(~) 


n =0 

oo 

(l- Z )-^J2g n (l-z) 4n/3 , 

n =0 
oo 

(l-z) 2 / 3 ^a n (l-z ) 4n/3 , 

n=0 


/o 00 

^ log (1-z) + g$ n (l-^) 4n /3 • 


(28a) 

(28b) 

(28c) 

(28d) 


We find that the entire function f(z ) is free in the IR, along with the values of B and <f>o- For 
simplicity, we fix f(z) = (1 — z) 2 / z 2 and choose B, tuning the value of <I?o so that the other metric 
functions give the correct AdS 4 behavior. In this case, we do not find the need to introduce a 
perturbation to allow us to flow to different values of the UV parameter that characterize the theory. 
Since here we do not have any electric component to our solution, we will take the dimensionless 
ratio to be B/<f> 2 . 

Figure [7(a) | shows the finite part of the entanglement entropy for the strip geometry for three 
solutions with different values of B / cj) 2 . We note that the behavior seems similar to that of the 
PED, suggesting that the form of the entanglement entropy may be invariant under electromagnetic 
duality. We explore this next. 



Figure 7: (a) The finite part of the entanglement entropy ss as a function of the strip width £ 
for the PMD with B/c /) 2 = 0.4183 (blue/lower), B/q if = 0.6710 (red/middle), and B / 4>\ = 1.6383 
(orange/upper), (b) The finite part of sd/£ as a function of l for the same value of the UV 
parameters with the same colors/locations. See footnote [6] for an explanation of the numerical 
variance in (b). 
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For the disk geometry, the finite part of the entanglement entropy can be found in Fig¬ 
ure |7 (b) | for the same values of the UV parameter. Again, the behavior is similar as to that of the 
PED disk entanglement entropy. 


6 Electromagnetic Duality 


For both the PED and PMD solutions, the Maxwell equation of motion from the action in eqn. © 


is 


d u 


—GZ (<£>) F^j = 0 . 


(29) 


In vacuum, there is the well-known electromagnetic duality where the equation of motion, eqn. (29), 
is invariant under F, IU —y *F, where *F is the Hodge dual given by 


(*r) u „ = S V=5, 


n px o auj F\ 

fiupcr^ ^ 1 X uj- 


(30) 


Here, £^ vpa is the completely antisymmetric tensor with £0123 = +1. In the case when there is 
a scalar in the theory that couples to the Maxwell field strength as in the action in eqn. © as 
Z(<&)F 2 , then the electromagnetic duality also requires that $( 2 ;) —> We now consider the 

electromagnetic duality in the boundary theory, applied to the PED and PMD solutions and their 
strip entanglement entropy. 


Evaluating eqn. (30) for the PED solution with metric written as in eqn. (21) we have that 

A(z) 


(*n* = - 


vn*)G(z 


h'{z) 


(31) 


Recall we require the UV behavior of h(z ) = /t — Qz and {F 1 , G, A} —> {1,1,1}. Thus, the UV 
“magnetic dual” solution to the PED solution with given Q value is 


(*F) xy = Q 


(32) 


Because we have scaled out the dependence of the solution on the horizon, zh, we must always con¬ 
sider dimensionless ratios. Thus in comparing the PED and PMD solutions under electromagnetic 
duality, we fix the value of Q/4>\ in the PED and then find the corresponding PMD solution that 
has B/(f)\ = Q/4> i- We also check that the UV behavior of the dilaton has opposite sign between 
the two theories. 

It turns out that all of our PED solutions have approximately the same value of Q ~ 1.224, 
although the dimensionless ratios change. We pick the value Q/4>\ = 3.40684, which is for a Sf 
value of 5f = —4.7218. We then find the corresponding PMD solution with B/(p\ = 3.40688. We 
note that for the PED solution, <j) 1 = 0.5995 and for the PMD </>i = —0.5994. We then compute 


the finite part of the strip entanglement entropy for both and obtain the results in Figure 8(a 
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(a) (b) 


Figure 8: (a) The finite part of the strip entanglement entropy ss(t) for the PMD (red/upper) and 
the PED (blue/lower), EM dual solutions with ratios B/cj>\ = Q/4>\ = 3.4068. The PMD is shifted 
up by 0.3440. (b) The same entanglement entropies as in (a) but with the PMD shifted down; see 
text for explanation. The agreement is nearly perfect. 


We find that the two entanglement entropies are shifted relative to one another, with the 
PMD solution being 0.3440 higher than the PED. We observed a similar shift between the analytic 
and numerical backgrounds for the PED solution Sf = 0, which was due to different coordinates 
giving rise to different UV cutoffs. For the PED and PMD solutions, we have used the same coor¬ 
dinates with the horizon at z = 1; however, there is a subtlety. We rescaled the radial coordinate z 
by the position of the horizon, zjj, so that we could use the range 0 < z < 1 for our numerics. 
However, in the definition of the entanglement entropy in eqn. this rescaling will change the 
UV cutoff e by e/ zh- Thus, if we are to restore the position of the horizon, we would have 


epED = £zq , (33a) 

epMD = eze , (33b) 


where e = 10 -5 is the numerical UV cutoff we have used in both backgrounds. If our choices of the 
PED horizon zq and the PMD horizon zb are not the same, it would lead to slightly different UV 
cutoffs between the two theories, and hence a shift in the entanglement entropy. 

To see this, recall from the analytic vs. numerical shift we observed previously, we know 
that the shift, rj, is related to the UV cutoffs via 


1 1 
= r] + 


epMD 


epED 


(34) 


Using (33), this means that the PMD horizon zb in terms of the shift ij and the PED horizon zq 
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is given by 


(35) 


ZB = 


Z Q 


1 + «7«Q 

Thus, if the shift is non-zero, then we have that zb ^ zq. We check this by shifting the PMD 
entanglement entropy down by ss(£) — 0.3440 for all values of £ and find near perfect agreement 
with the PED entanglement entropy - see Figure [8(b)[ 

As a further check, we revisit the analytical solution m for the PED, given in our notation 
by eqns. (24). We can now do a rescaling of the analytic solution’s coordinates, letting z' = z/Q, 
which mimics the rescaling of our coordinates by the horizon radius. This implies that the UV 
cutoff e' for the analytic solution is then rescaled to become e/Q and we find a relationship similar 
to eqn. (33): e = eQ. Following the above arguments for our case with the different horizons, 
we find that if Q is varied then there will also be a shift r / in the analytic entanglement entropy, 
controlled by Q, via 


Q = 


l 


1 + e'rf 


(36) 


This shows that changing Q should have the same effect as changing the horizon positions in 
the PED and PMD solutions, which leads to a shift in the entanglement entropy. We show in 

entanglement entropy ss{£) for values of Q = 1,1.5, 2 and indeed there 
we subtract this shift from the entanglement entropies when Q > 1 and 
we find there is good agreement, which confirms that changing Q-ancl hence changing the horizon 
position for the PED and PMD solutions-leads to an overall shift. 


Figure 9(a) the analytic 


is a shift. In Figure 9(b) 


SS ss 




Figure 9: (a) The finite part of the strip entanglement entropy ss(£) for the analytic PED solution. 
The three curves are for Q = 1 (blue/lower), Q = 1.5 (red/middle), and Q = 2 (orange/upper), 
(b) The same entanglement entropies as in (a) after being shifted down; see text for explanation. 

From the point of view of the dual field theory, the electromagnetic duality should act 
by exchanging electric and magnetic charge carriers [371 EH], and so our results suggest that the 
entanglement entropy is invariant under this exchange, at least for the field theories dual to the 
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PED and PMD background^} 


7 Entanglement Entropy of the Dilaton—Dyon Black Hole 


Now we turn to the dilaton-dyon black hole. 

The IR series expansion for the metric in 

eqn. (fTlj) 

Maxwell field (12 

), and dilaton are 





F(z) 

= (1- 

OO 

n= 0 

(37a) 


G ( z ) 

1 

" (1- 

OO 

n—0 

(37b) 


M*) 

iWl 

II 

n(l ~z) n , 

(37c) 


Hz) 

— ho(l 

~ z) , 

(37d) 


Hz) 

= I>n(l ~z) n . 

(37e) 


71=0 


The equations of motion fix all coefficients except for {/o,ai,3>o}> and we also have a choice of 
the value of the magnetic held B. We fix /o = 1 and generate different backgrounds that give 
different dual UV held theory parameters by choices of B,a\ and $o- Our particular choice of 
ansatz actually hxes Q = B for all of our solutions. 

In Figure 10(a), we plot ss for the dilaton-dyon black hole for a few values of B / fj?. We hnd 
that increasing the value of B/fi 2 appears to shift the entire entanglement entropy down. Zooming 
in on ss(£) (Figure 10(b)), we see that the behavior of the entanglement entropy appears to be 
linear in £, ss(£) ~ £, with the slope increasing as B / fi 2 is decreased. 

The results for the dislj^] entanglement entropy are shown in Figure 11, where we again 
found it useful to switch to the radial coordinate dehned in (27), for three values of B/fi 2 . We 
hnd that the entanglement entropy, Sr> /^ , appears to be linear in £, which suggest that sd{£) ~ £ 2 , 
which would be a “volume” law rather than an area law for the entanglement entropy. As with 
the strip entanglement entropy, increasing B / fi 2 shifts sr>/£ down but it also leads to a much more 
noticeable shift in the slope of the linear part of sd /F 

As a check, we turn off the dilaton (consistently setting <h(z) = 0 in our equations of motion) 
and compute sd/£ for the dyonic black hole. We can compare our numerical background’s result 


'In ref. [27], they also find that under electromagnetic duality transformation, the logarithmic violation of the 
area law that they see is preserved, suggesting a Fermi surface even though the charge density vanishes. 


8 


The numerical results for the dilaton-dyon sd in Figure 11 have very little numerical variance in comparison to 


those for the PED and PMD solutions. This is due to the fact that for the dyon solution the dilaton does not diverge 
in the IR. 
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Figure 10: (a) The finite part of the entanglement entropy ss as a function of the strip width t for 
three values of the dimensionless ratio B/^i 2 . Blue/upper curve is B / [i 1 = 1.040, red/middle curve 
is B / jj? = 6.779 and orange/lower curve is B / fi 2 = 27.866. (b) The same plot as (a) but zoomed 
in to see the large l behavior of ss- 




Figure 11: (a) The finite part of the entanglement entropy sd/1 as a function of the disk radius 
l for four values of the dimensionless ratio B/fi 2 . Blue/upper curve is B / fi 2 = 3.249, red/middle 
curve is Bf\i 2 = 6.779, and orange/lower curve is £>//i 2 = 27.866. (b) The same data as in (a) but 
plotting sd/I 2 . See footnote [8] for an explanation as to why the results have less numerical variance 
than those for the PED and PMD. 


to that of the analytical background written in ref. [39]. The result is shown in Figure 
see that there is very good agreement. 

All of our solutions are at zero temperature, and the PED black holes also have zero thermal 
entropy because they have zero area horizon, but our dilaton-dyon black hole has a finite area 
horizon (which in the dual gauge theory is interpreted as a large ground state degeneracy) and thus 
a finite entropy, S, even at zero temperature. Namely, computing the area of the horizon at z = 1 


12 and we 
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Figure 12: Comparing the numerical and analytical dyonic black hole sd/£ 2 . Blue is the analytical 
background and red is the numerical background. 


Plot Color/Location 

UV Bln 2 

IR B 

B 1 / 2 / (3 1 / 4 • 4) 

Slope 

Blue/Upper 

3.249 

0.3074 

0.1053 

0.1011 

Red/Middle 

6.779 

0.1477 

0.0730 

0.0569 

Orange/Lower 

27.866 

0.035 

0.0359 

0.0322 


Table 1: Comparing the numerically determined slope of the linear part of Figure [1 1(a) for sr>/t 
to the ground state degeneracy in eqn. (39). The plot colors/location and UV value of B / fi 2 


are 


given for ease of reference, as is the IR value of B , which is used to determine the ground state 
degeneracy via eqn. (39). 


using our IR expansion (37) (and dividing out by the infinite volume of M 2 ), we find 


S 

Vol(M 2 ) 


Area(Horizon) 

4 



(38) 


which is a constant, where we have absorbed the factor of the AdS radius L into the volume factor. 
In fact, from the equations of motion, we have that qq = and so 


S _ 1 B 1 / 2 
Vol(K 2 ) “ 3V4^T 


(39) 


To see whether or not our entanglement entropy approaches the ground state entropy in 


limit, we did a numerical fit to the linear part of sjj/£ in Figure 11(a) The resulting 
Table [IJ comparing to the value of eqn. (39) using the IR value of B. 


the large i 
slope is in 
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8 Entanglement Entropy of the Magnetic Electron Star 


We are now ready to turn on the star. The solutions were found in ref. [26] in coordinates where 
the horizon is at infinity, but here we find it better to rescale those solutions so that the horizon is 
at z = 1. We will consider the mesonic phase, where there is a charged star in the IR as well as a 
horizon, but with no electric charge behind the horizon. The IR expansion of the fields is given by 


F(z) = z 2 (l- z) 2 f n { 1 - z) 2n/3 + <5/(1 - z ) 2b/3 , 

\n =0 ) 

G{z) = z 2 { 1 - z)-^ (x>n(l - ^) 2n/3 + <5*7(1 - *) 26/3 ) , 

\n =0 ) 

A{z) = z 2 ( 1 - z ) 2 / 3 ( ^a„( 1 - z ) 2n / 3 + <5a(l - z) 2b ^ J , 


\n=0 


h(z) = h 0 ( 1 - z), 

/q 00 

$(z) = ^ In (1 -z) + J2 $n(l - z) 2n/ 3 + < 5$(1 - z) 2b / 3 . 


(40a) 

(40b) 

(40c) 

(40d) 

(40e) 


n=0 


We have again turned on a perturbation with b = | (—3 + \/57) in order to flow to different values 


of the UV parameters (14). The equations of motion fix all coefficients except {!?,/o, •ho, h*!*}. We 


choose a value of the fermion mass, rh, and choose { B , /o, and then adjust 4>o (and rescale the 


time coordinate) so that the solution flows to AdS 4 in the UV. Figure 13(a) shows an example, 
plotting the local chemical potential, fermion mass, and value of the ending radius of the star, z s t a r\ 


Figure 13(b) shows the metric functions and dilaton for the same background. 


For this background, we show in Figure 14(a) the finite entanglement entropy for the strip 


geometry. As with the PED black hole solution, we have computed the area of the “in-falling” 
solution. We find that the large £ value of sg is greater than the area of the “in-falling” solution. In 
fact, we are able to generate several magnetic electron star backgrounds (see Figure [15]) and from 
those backgrounds there are strong indications that for every one the large l value of sg is greater 
than the in-falling solution, with the difference between the two always around 1.5 ±0.2. However, 
we do not show curves for sg for all of these backgrounds because the data are too difficult to 
control numerically well enough to display - there is too much numerical noise; nevertheless, the 
relationship between the in-falling and large £ value of sg is clear to us from these investigations. 
We also checked this behavior in a different coordinate system (one in which the horizon is at 
z —> oo) and again find the same behavior. This suggests that the behavior is robust. Since 
the in-falling solution has a lower entropy than the minimal surface, we conclude that the actual 
entanglement entropy as seen in Figure [14] follows the red curve, then the blue, until it meets 
the orange line, and then it remains that of the in-falling solution for all £ thereafter (that is, 
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Figure 13: (a) The local chemical potential is shown in blue (the curve), where the fermion mass 
(shown in red, a horizontal line) is rh = 0.5 and z s t ar = 0.3529 is the vertical line. The dual UV 
field theory parameters are B/g 2 = 6.8132, Q/ /i 2 = 0.1136, 4>i/g = —3.6507, and fo/£? = 12.4852. 
(b) The background metric and dilaton (lowest purple curve) for the same UV parameters as in 
(a). Here we have actually plotted G(z) = z 2 g(z), A(z) = z 2 a(z), and F(z) = z 2 f(z) (lower-middle 
blue curve, upper red curve, and upper-middle orange curve respectively) so that our expect UV 
behavior is that each metric function goes to 1. We have also pulled out a factor of 2 in ^(z), 
redefining <f>(,z) -» z$>(z) (purple). The vertical line is z s t ar - 


SS SS 



Figure 14: (a) The finite part of the entanglement entropy ss for the strip geometry as a function 
of the strip width t for the mesonic phase of the magnetic electron star with B/g 2 = 6.8132. The 
red, lower part of the curve (that tends to negative infinity), is the part of ss that lies outside the 
star and the blue, upper part of the curve (larger values of £), is the part that lies inside the star; 
the minimal surface crosses into the star at a value of £ s t ar = 0.1938. The orange, horizontal line 
is the value of ss for the “in-falling” solution and the intersection of the horizontal and vertical 
lines indicate the point at which the entanglement entropy crosses over to the in-falling one. (b) 
Same as (a) but for B/g 2 = 14.1343. Here the minimal surface crosses into the star at a value of 
^star = 0.1883. This differs from (a) since the value of z s tar is different. 
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it follows the curve starting from negative infinity all the way until it means the intersection of 
the vertical and horizontal lines, after which it remains the horizontal line). From the point of 
view of the minimal surface, at a certain value of i, the hanging surface breaks and becomes the 
in-falling solution. Figure [l4(b)| shows the strip entanglement entropy for a different value of the 
dimensionless ratio B/^j, 2 to demonstrate this point. We also find evidence that for the backgrounds 
generated along the line in Figure 15 that the entanglement entropy increases as B / fj 2 is increased 
(as explained above, we choose not to display this data due to the difficult numerics involved). 


Q_ 

a 2 



ti 

A 



B 

A 2 


Figure 15: (a) The line in the phase diagram for the different magnetic electron star backgrounds 
we are able to generate, showing here the relationship between B/fi 2 vs. Q//r 2 , where B/[i 2 is the 
free parameter, (b) The same line in the phase diagram but showing two of the parameters we 
control, vs. B/[i 2 . 


We note that we were unable to compute the disk entanglement entropy for the mesonic 
phase of the magnetic electron star at this point in time due to numerical instability in solving the 
second order problem for these backgrounds, at least at the level of the numerical sophistication 
that we are using. In addition, although we have solutions for partially fractionalized phases of 
magnetic electron stars, given in ref. [26], we are currently unable to change how populated the 
star is (i.e. how much the chemical potential is greater than the fermion mass) and so cannot fully 
explore the entanglement entropy of these solutions. 

Our results strongly indicate that the entanglement entropy of the mesonic phase of the 
magnetic electron star is undergoing a “phase transition” as it transitions to the in- falling solution. 
This is reminiscent of the breaking of the hanging string in confinement/deconfinement phase 
transitions mi ehe], and similar behavior has been seen in the work of two of the authors of 
this paper, ref. [323’ hi the case of holographic superconductors. However, without being able to 
compare to the disk case we cannot say how robust this “phase transition” may be, and in addition, 
without a careful study of the stability of magnetic electron stars we cannot rule out the possibility 
that we are detecting an instability. These are matters that demand further investigation. 
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9 Conclusion 


We have computed the finite part of the holographic entanglement entropy for the strip and disk 
geometry for purely electric, purely magnetic, and dyonic dilatonic black holes, as well as the strip 
entanglement entropy for the magnetic electron star solutions of ref. [26]. We observed how the 
entanglement entropy for the PED solution backgrounds at large l approach the in-falling solution, 
which contributes a non-zero value to the entanglement entropy. We also found that different UV 
regulators, due to having different coordinate systems, can introduce a shift in the entanglement 
entropy. Similar results are obtained for the PMD backgrounds, and we are able to show that the 
entanglement entropy is invariant under electromagnetic duality, as long as one shifts the relative 
location of the physical horizons between the PED and PMD solutions. 

For the dilaton-dyon black hole, we found that the entanglement entropy for the disk, at 
large £, grows linearly with £, with a coefficient approximately equal to the ground state degeneracy 
of our solutions. This suggests a volume law rather than an area law for these backgrounds. Lastly, 
for the mesonic phase of the magnetic electron star solution, we found that the entanglement 
entropy smoothly passes into the star, but that at a certain value of £, it is actually the in-falling 
solution that has lower entropy than the hanging surface. This suggests that at a value of z, with 
z s tar < z < 1, the minimal surface breaks to the in-falling one. This strongly suggests the presence 
of a “phase transition,” however, pending further investigation we cannot say for certain; regardless, 
our results indicate at least some non-trivial behavior in the far IR for the entanglement entropy. 
In all cases, we find a dependence of the large £ behavior of the entanglement entropy on the various 
dimensionless ratios in the theory. 

There are several directions for future work. We mentioned in Section [8] that one could 
study the disk entanglement entropy for the mesonic case, as well as the entanglement entropy of 
the partially fractionalized phases of the magnetic electron stars. For the case of zero magnetic 
field but with a charged star and electrically charged horizon-a partially fractionalized phase-the 
entanglement entropy’s dependence on the charge was studied in refs. [12] and there it was found 
that there was a logarithmic violation to the area law due to the charged horizon but not the 
star. It would be interesting to see if the presence of the magnetic held in such phases changes 
this behavior. Since we have only studied the mesonic phase in this paper, we do not expect 
the logarithmic violation, at least in the electric sector. Perhaps, however, even in the mesonic 
phase the presence of the horizon for the magnetic sector will modify the area law behavior of the 
entanglement entropy in an interesting way. Additionally, one could construct and study finite 
temperature versions of our solutions. These points are left for future work. 
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Appendix 

We present here the details of the derivations of eqns. <©> 0, and ([8]) given in Section [2| Recall 
we have taken bulk metric to be of the form 


ds 2 = -g a (z)dt 2 + g zz (z)dz 2 + g xx (z)dx 2 


(41) 


where UV is at z = 0 and the horizon at z = 1, after rescaling the coordinates by zh, the physical 
horizon position. Writing the M 2 part of the metric as dx 2 + dy 2 and having the finite length of the 
strip run from —1/ 2 < x < t/ 2, the minimal surface will be symmetric about y = 0. We denote 
the z value of the turning point of the surface as zt- Let the surface have coordinates (x,y), the 
same as the M 2 coordinates, with embedding z = z(x). Then we may write the pull-back of the 
metric, at fixed time, as 


ds iA = ( 9zz{z{x)) ( ^ 


+ 9xx(z(x)) dx 2 + g xx (z(x))dy 2 


(42) 


and thus the finite part of the entanglement entropy becomes 


4G N S s = 2Jz? 


r t /2 


dxg xx (z(x )) 



9zz{z{x)) ((fey 
9xx(z(x)) \dxj 


(43) 


Because the above equation has no explicit dependence on the coordinate x, we can reduce the 
second order problem of minimizing this integral for the surface z(x) into a first order problem M- 


We see that eqn. (43) does not have any explicit dependence on x, so thinking of it as a time 
coordinate the associated “Hamiltonian” will be a conserved quantity with respect to x. i.e., = 0. 


We denote by p z the conjugate momentum for the variable z, and taking the integrand in eqn. (43) 
as the “Lagrangian”, we find 


Pz = 


9zz(z(x))z'(x) 


1 + 


9zz(z(x)) 


gxx{z{x)) 

which allows us to write the conserved “Hamiltonian' 


(z'(x)Y 


(44) 


H = - 


gxx(z(x)) 


1 + 


9zz{z(x)) 

gxx(z(x)) 


= constant. 


(z'(x)y 


(45) 
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This constant of motion is related to the turning point zt, as we will shortly see, so we denote the 
constant by c(zt) and, squaring both sides of eqn. (45) and doing some algebra, we find 


dz 

dx 


1 Qxxiz') 
c 2 (z t ) g zz {z) 


(9%x(z) - c2 ("T)) 


(46) 


This is our first integral of motion, reducing the original second order problem of finding z(x) to 
this first order one. We can go further, however, and bypass solving for z(x) directly. First, from 


eqn. (46) we see that the constant c{zt) is related to the turning point if we set 


c{z T ) = ±g xx (z T ) • 


Then we have that 


dx = 


g xx (z T )dz 


9xx (z) 


9zz 


(ij (9xx( z ) - 9 xx(zt)) 


(47) 


(48) 


This last relationship allows us to write the original entanglement entropy integral eqn. (43) as an 
integral over z. The result is 


4 G n S s = 2«£f 


fZ T 


dz 


g xx {z) 


\J Em ( 9x X {z ) 2 - g X x(z T ) 2 ) 


(49) 


where we have introduced the UV cutoff e in the lower limit of the integral. With this, we can 
easily find the length of the strip, t = 2 f^ 2 dx, as a function of the turning point via substitution 
of dx above; we find 

g xx (zt ) 


rz T 


' dz — = . 

0 vEm (sxx(z) 2 - g xx (z T ) 2 ) 


(50) 


Since the length t should be positive, in writing eqn. (50) we made the choice c(zt ) = —g xx {zT)- 


The same trick does not work, however, for the disk geometry, and so the formula is simply 
given by the pull-back 


4:G n S d = 2it[ drrg xx (z(r))\ 1 + 9zz ^^ (—\ 


(51) 


g xx {z{r)) \dr J 

where 0 < r < t is the radial variable for polar coordinates of the M 2 , dr 2 + r 2 d0 2 . In this case, 


we must explicitly solve for the minimal surface z(r) and then input that into eqn. (51) to find the 
entanglement entropy. Eqns. (49), (50), and (51) are the same as eqns. ([6]), ([7]), and ^ given in 
Section [2j 
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